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We propose an approximation to the ground state of Yang-Mills theory, quantized in temporal gauge and 2+1 
dimensions, which satisfies the Yang-Mills Schrodinger equation in both the free-field limit, and in a strong- 
field zero mode limit. Our proposal contains a single parameter with dimensions of mass; confinement via 
dimensional reduction is obtained if this parameter is non-zero, and a non-zero value appears to be energetically 
preferred. A method for numerical simulation of this vacuum state is developed. It is shown that if the mass 
parameter is fixed from the known string tension in 2-1-1 dimensions, the resulting mass gap deduced from the 
vacuum state agrees, to within a few percent, with known results for the mass gap obtained by standard lattice 
Monte Carlo methods. 



I. INTRODUCTION 

Confinement is a property of the vacuum state of quan- 
tized non-abelian gauge theories, and it seems reasonable that 
something could be learned about the origin of confinement, 
and the origin of the mass gap, if we knew the form of the 
Yang-Mills vacuum wavefunctional in some physical gauge. 
There have, in fact, been a number of efforts along those lines, 
in temporal gauge [1-7], Coulomb gauge [8, 9], axial gauge 
[10], and in a Bars corner-variable formulation [11, 12]. 

In this article we will pursue this investigation in tempo- 
ral gauge and in D = 2 + 1 dimensions, our strongest influ- 
ences being refs. [1] and [7]. Our claim is that the ground 
state wavefunctional ^o[A] can be approximated by the form 



^o[A] = exp 



-^J cfx(fy B"{x) 



ab 



B'iy) 



(1) 



where B" — ¥"2 is the color magnetic field strength, — 
Dj^Dii is the two-dimensional covariant Laplacian in the ad- 
joint color representation, A-o is the lowest eigenvalue of —D^, 
and m is a constant, with dimensions of mass, proportional to 
g^. To support this claim, we will argue that the above expres- 
sion 

1. is the ground state solution of the Yang-Mills 
Schrodinger equation in the g — > limit; 

2. solves the zero-mode Yang-Mills Schrodinger equation 
in the zero-mode strong-field limit; 

3. confines if m > 0, and that m > is energetically pre- 
ferred; 

4. results in the numerically correct relationship between 
the mass gap and string tension. 

A very similar proposal for the vacuum wavefunctional, with 
A-o absent, was put forward by Samuel in ref. [7], generalizing 
the earlier "dimensional reduction" proposal of ref. [1]. 



Our paper is organized as follows: In section II, below, 
we find an approximate solution of the zero mode Yang-Mills 
Schrodinger equation in 2h-1 dimensions, and compare this to 
our proposed wavefunctional in an appropriate limit. The di- 
mensional reduction and confinement properties are discussed 
in section III. Section IV outlines a procedure for numerical 
simulation of our vacuum wavefunctional; in section V this 
procedure is applied to calculate the mass gap, with param- 
eter m chosen to give the coiTect string tension as a function 
of coupling. Confinement, in our approach, relies on > 0; 
in section VI we will discuss why this choice lowers the vac- 
uum energy in the non-abelian theory, while the minimum is 
at = in the free abelian theory. Section VII contains a 
few results and critical comments regarding certain other pro- 
posals for the Yang-Mills vacuum wavefunctional. Some brief 
remarks about Casimir scaling and N-ality are found in section 
VIII, with conclusions in section IX. 

We would like to note here that the work in section II, con- 
cerning the zero-mode strong-field limit, was motivated by a 
private communication from D. Diakonov to one of the au- 
thors [13]. 



II. THE FREE FIELD AND ZERO MODE LIMITS 

In temporal gauge and D — d+l dimensions, the problem 
is to find the ground state of the Yang-Mills Schrodinger equa- 
tion 



where 



52 



(2) 



(3) 



and all states in temporal gauge, in SU(2) gauge theory, are 
subject to the physical state condition 



(4) 



This condition requires invariance of ^[A] under infinitesimal 
gauge transformations. 
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Our proposed vacuum wavefunctional, eq. (1), obviously 
satisfies the physical state condition, since the kernel 



^xy — 



1 



(5) 



xy 



transforms bilinearly, Kxy U{x)KxyU~^{y), under a gauge 
transformation, with U a transformation matrix in the adjoint 
representation. In the g — > limit, with both A,o,ot — > in the 
same Umit, the vacuum state becomes 



exp 



sab 



dhd^y(diA^ix)-d2A"i{x)^ 
(d,A'2{y)-d2A\{y)) (6) 



which is the known ground state solution in 2+1 dimensions, 
in the abelian, free-field case. 

The Yang-Mills Schrodinger equation is also tractable in a 
quite different hmit, which is, in a sense, diametrically op- 
posed to the free-field situation. Let us restrict our attention 
to gauge fields which are constant in the two space directions, 
and vary only in time (analogous to the minisuperspace ap- 
proximation in quantum gravity). The Lagrangian is 



/ 



d,Ak ■ dtAk - g^{Ai X A2) • (Ai x A2) 



2^ 



dtAk ■ dfAk - g^{Ai X A2) • (Ai x A2) 



(7) 



it is not hard to verify that 



To = 0+-g^V 



IVRl 



(Ai XA2) • (Ai x A2) 



(|AiP + |A2|2)' 



|Ai|2 + |A2p 



(14) 



Now for A-fields for which *Po is non-neghgible, it is easy to 
see that TqI'o is of order no greater than 1 /V, except in the 
immediate neighborhood of the origin (A^. = 0) of field space. 
That is because ^0 ~ exp[— y/?o], which is non-negligible 
only if y/?o is 0{\). For comparison with eq. (1) we are inter- 
ested in a strong-field limit, far from the origin of field space. 
In that case, since /?o ~ 1 /V, then the rhs of (14) is at most of 
order 1 /V, which can be neglected. It follows that Rq in eq. 
(12) accomphshes the required cancellation at leading order, 
and provides the leading contribution to the logarithm of the 
vacuum wavefunction. 

Now consider the proposal (1) for the vacuum wavefunc- 
tional of the full theory, in a comer of field space where the 
non-zero momentum modes of the A-field are negligible com- 
pared to the zero modes, and in fact the zero modes are so 
large in magnitude that we can approximate D'j^ w ge'^^'^A*. 
In this region 



{-DX = g'?>\x-y)M' 



■ab 



(15) 



where 



where V is the area of a time-slice, leading to the Hamiltonian 
operator 



H 



The factors of V in the Hamiltonian suggest the use of a 
1 /V expansion. Let us write 



n = exp[-y/?o + V"^/?2 + ...) 



(9) 



with Ro chosen such that the leading order (in 1 /V) "kinetic" 
term contained in //^o 
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cancels the potential term 
atO(y). Let 



i^2y(AixA2)-(AixA2)1'o 



J (Ai XA2) • (Ai X A2) 
«o = 2S 



Then, defining 
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In SU(2) gauge theory, the two zero-mode fields Ai , A2 define 
a plane in three-dimensional color space. Take this to be, e.g., 
the color x — y plane, i.e. 



(17) 





a\ 




bi 


Ai = 


fl2 


, A2 = 


bi 













Then 



M = 



aj + bl 
-aia2 — bib2 




-(31(32 — bib2 
a\ + b\ 






A\+Al 



Now M has three eigenstates 
^1 











"0" 




, ^2 = 




, <l>2 = 















1 



with corresponding eigenvalues 



A'2 



(18) 



(19) 



(20) 
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where 



Then 



S=aI+aI , C={AixA2)-{AixA2) 



ab 
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(21) 
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We have 



exp 




= exp 

{AixA2f 

Taking account of eqs. (17), (19) and (22), we get 

1 



(23) 



exp 
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-\gV{AixA2f 



(AixAz)^ 



= exp 



1 (Ai xA2)-(Ai XA2) 



Now by assumption, in the strong-field hmit. 



g'Hi^g\A\+Ai):^m' 



(24) 



(25) 



and 



W = 5^ l-\/l-4 



C 
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We recall that the ground- state solution of the zero-mode 
Schrodinger equation ^0 ~ exp[— y7?o] with Rq given in eq. 
(12) is vahd for Rq ~ l/V, where the wavefunction is non- 
negUgible. In this same region of configuration space, /ii is 
negUgible compared to /i3, and eq. (24) becomes 



exp 



I (Ai X A2) ■ (Ai XA2) 
' 2^^ 



VM3 



= exp 



j_ (Ai xA2)- (Ai XA2) 



(27) 



which is identical to the solution found for the ground state 
of the zero-mode Schrodinger equation, in the region of va- 
lidity of that solution, where VRq ~ 0(1) . Therefore, in a 



small region of configuration space where a non-perturbative 
treatment is possible, we find that our ansatz for the vacuum 
state agrees with the ground state of the zero-mode Yang- 
Mills Schrodinger equation. ' 

The argument above can also be extended to 3+1 dimen- 
sions, as outhned in Appendix A. 



m. DIMENSIONAL REDUCTION AND CONFINEMENT 

Assuming that our proposal (1) for the Yang-Mills vacuum 
wavefunctional in 2+1 dimensions is at least approximately 
correct, then where does the confinement property appear? 

A long time ago it was suggested that the effective Yang- 
Mills vacuum wavefunctional at large scales, in D = d + 1 
dimensions, has the form [1] 



•' exp 



(28) 



(see also [3, 5]). This vacuum state has the property of "di- 
mensional reduction": Computation of a large spacelike loop 
in d + 1 dimensions reduces to the calculation of a large Wil- 

(3) 

son loop in d Euclidean dimensions. Suppose 'Pq ' is the 

(2) 

ground state of the 3+1 dimensional theory, and ' is the 
ground state of the 2+1 dimensional theory. If these ground 
states both have the dimensional reduction form, and W(C) 
is a large planar Wilson loop, then the area law falloff in 
Z) = 3 + 1 dimensions follows from confinement in two Eu- 
cUdean dimensions in two steps: 



W(C) = (Tr[f/(C)])«=4 = (<)|Tr[f/(C)]|^^^^) 



/(3)\ 



(Tr[f/(C)])°=3 = (^(2yr[[/(C)]|^^^^) 
(Tr[[/(C)])^=2 



/(2)\ 



(29) 



In Z) = 2 dimensions the Wilson loop can of course be cal- 
culated analytically, and we know there is an area-law falloff, 
with Casimir scaling of the string tensions. The dimensional 
reduction form of the ground state wavefunctional can be 
demonstrated explicitly in strong-coupling lattice gauge the- 
ory [2]; Monte Carlo support for the hypothesis has also been 
obtained at intermediate couplings [14, 15]. 

It is natural to try and improve on the dimensional reduc- 
tion idea by considering wavefunctionals which interpolate, 
in some natural way, between free-field dynamics at short 
distance scales, and the dimensional reduction form at large 
scales. In ref. [7J, Samuel suggested that the vacuum state in 



' Wc learned from D. Diakonov that he had obtained this result in unpub- 
lished work, which considered a wavefunctional of similar form to (1) but 
without the ^o,'« terms in the kernel [13]. Those terms are not important 
in the region of configuration space discussed in this section. 
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D = 2+1 dimensions might have the form 



or in lattice units, with lattice coupling p, 



'Po[A]=exp 



ab 



£)2 



xy 



(30) 

Our proposal differs from Samuel's in that wiq is replaced 
by — Xo + m}, with the lowest eigenvalue A,o being field- 
dependent and gauge-invariant. The rationale is that we 
should allow for a subtraction in the operator appear- 
ing in the vacuum kernel; a subtraction will be absolutely re- 
quired if the spectrum of —D^, starting with Xq, diverges in 
the continuum limit. On the other hand, if < is a nega- 
tive constant, then the wavefunctional in eq. (30) is not neces- 
sarily real throughout configuration space, and can oscillate. 
Now the true vacuum state must be real up to a constant fac- 
tor, and it is forbidden to pass through zero by the "no node" 
theorem for quantum-mechanical ground states. Requiring a 
subtraction which respects the reahty of the wavefunctional, 
and avoids oscillations anywhere in field configuration space, 
dictates the replacement 



ml- 



-'ko + m 



(31) 



with > 0. 

The dimensional reduction form is obtained by dividing the 
field strength into "fast" and "slow" components, defined in 
terms of a mode cutoff. Let {(|)"} and {%„} denote the eigen- 
modes and eigenvalues, respectively, of the co variant Lapla- 
cian operator in adjoint color representation, i.e. 



The field strength can be expanded as a mode sum 

B''{x)=j^bnm 
n=0 

and we define the "slow" component to be 

ihnax 

n=0 



(32) 



(33) 



(34) 



where rimax is a mode cutoff chosen such that AX = X„„^ — 
A-o <C in^ remains fixed as V ^ °°. In that case, the portion 
of the (squared) vacuum wavefunctional gaussian in is 
approximately 



exp 



L /" Qa,slowQa,slow 

m J 



(35) 



which is just the probability measure for Yang-Mills theory in 
two Euclidean dimensions, with a particular type of ultravio- 
let cutoff. The string tension for fundamental representation 
Wilson loops in D = 2 Yang-Mills theory, with couphng g^m, 
is easily computed: 



o=-8m 



(36) 



3 m 



(37) 



In the next sections we will address two questions. First, 
suppose we fix m to give the known string tension at a given 
lattice coupling. What is then the value of the mass gap pre- 
dicted by the vacuum wavefunctional, and to what extent does 
this agree with the corresponding value determined by stan- 
dard lattice Monte Carlo methods? Secondly, since confine- 
ment depends on having w 7^ 0, is there any reason why the 
mass parameter m should be non-zero? 



IV. NUMERICAL SIMULATION OF THE VACUUM 
WAVEFUNCTIONAL 

The mass gap implied by the vacuum state (1) can, in prin- 
ciple, be extracted from the equal-times connected correlator 

©(x-y) = {{B"ir'UB''B'')y) - {{B-B"),)^ (38) 

where the expectation value is taken with respect to the proba- 
bility distribution P[A] defined by the vacuum wavefunctional, 
i.e. 

(2) = J DA1DA2 Q[A]P[A] (39) 



with 



P[A] = \^o[A]\ 



= exp 



~ f d^xcfyB%x)K-^\A]B'^{y) 



and 



ab 



(40) 



(41) 



Here we have absorbed a factor of g into the definition of A,, 
which accounts for the factor of l/g^ in the exponent in eq. 
(40). 

It not easy to see how ® (x — y) could be computed analyt- 
ically beyond the level of weak-coupling perturbation theory, 
but computation by numerical simulation of P[A] also seems 
hopeless, at least at first sight. Not only is the kernel 
non-local, it is not even known explicitly for arbitrary A^(x). 
However, suppose that after eliminating the wild variations of 
K along gauge orbits via a gauge choice, K[A] has very lit- 
tle variance among thermalized configurations. In that case, 
things are more promising. 

Let us define a probability distribution for gauge fields A 
which is controlled by a second, independent configuration A' 



A;K[A' 
detV2 



exp 



(42) 
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where the field strength B is computed from the A- 
configuration, and both A and A' are fixed to some appropriate 
gauge. Now, assuming that the variance of K[A] in the proba- 
bihty distribution P[A\ is small after the gauge choice, we can 
approximate 

P[A\ « p[a,{K)\ 

= p]^, j DA' k[a!]p[a!] 

« j DA' p'^,K[A'] P[A!] 



(43) 



where the step from the second to the third line follows from 
assuming that the variance of K in the distribution P[A\ is 
small. If this assumption about K[A\ is correct, and eq. (43) 
holds, then the probability distribution could in principle be 
generated by solving (43) iteratively: 



/'(i)[A] 



A;K[Q\ 



p(''+i)[A] = j DA' p\a-K[A']^P^''\a!] (44) 

A numerical version of this approach would be to use 
equibiUbrium configurations of /"("'[A], generated at the n-th 
step, to generate equiUbrium configurations of at the 

(«+l)-th step. 

We may use the remaining gauge freedom to fix to an axial 
gauge in the two-dimensional time slice. This allows us to 
change variables in the functional integral over two-dimension 
configurations from Aj, Aj to B^, without introducing a field- 
dependent Jacobian. Let eigenvalues and eigenmodes 
solve the eigenvalue equation 



(45) 



for the covariant Laplacian —D^ determined from the fixed A' 
configuration, and let {/7„} be the mode amplitudes of the B- 
field, as seen in the mode expansion (33). Then the probability 
distribution for the {bn}, which follows from P[A;A'[A']] at 
fixed A', is Gaussian 



prob[fc„] oc exp 



4 y'Xn — 'ka^ 



(46) 



In practice we use a lattice regularization on an L x L lat- 
tice with periodic boundary conditions, and the gauge field 
A^(x,>') is initialized to zero at the first iteration. We then 
generate gauge fields recursively; the procedure at the n-th it- 
eration is as follows: 

1. From one of the lattice configurations generated at the 
(n — l)-th iteration, compute the link variables in the 
adjoint representation, and then determine numerically 
the eigenvalues and eigenmodes of the two-dimensional 
lattice covariant Laplacian operator —D^. 

2. Generate a set of 3L^ normally-distributed random 
numbers with unit variance, denoted {r„}. From these. 



we obtain a new set of mode amplitudes 



and a corresponding B-field 

3Z.2 

B"{x) = Y,bnm 

n=0 



(47) 



(48) 



From the field strength fi"(x), and the axial gauge con- 
dition, determine the corresponding gauge field Af{x). 
This step can be repeated to generate as many thermal- 
ized configurations of /'[A,/r[A']] as desired. 

3 . The gauge fields are exponentiated to give link variables 

Uk{x,y)=exp[iAt{x,y)Ca/2] (49) 

and any observables of interest are computed. This con- 
cludes the n-th iteration. 

Lattice configurations generated by this procedure will be re- 
ferred to as ''recursion lattices". 

Details about our particular choice of axial gauge on a finite 
lattice, and the procedure for obtaining the A-field from the B- 
field in that gauge, may be found in Appendix B. 

V. THE MASS GAP 

The simulation procedure outhned in the last section leans 
heavily on the assumption that there is little variance in the 
kernel in a fixed gauge, or, equivalently, that there is neg- 
ligible variance, among thermalized configurations, in gauge- 
invariant combinations of the kernel such as Tr[A'~'^r~'], or 
in the gauge-invariant spectrum of K. The absence of signif- 
icant fluctuations in these quantities, when evaluated numeri- 
cally, is a self-consistency requirement of the method we have 
proposed. The quantity Tr[/Qy'/Cj;^'] is of particular interest, 
because its rate of falloff at large \x — y\ is determined by the 
mass gap. 

We begin with the spectrum {X„ — ^ -f w^} of the operator 
—D^ — T^+rrP', with m chosen, at a given P, to reproduce 
the string tension a(P) known from Monte Carlo simulations 
of the standard Wilson action in three Euclidean dimensions 
[16]. From eq. (37), this means choosing 



m = -Pa(P) 



(50) 



The result for the spectrum at (3 = 18 on a 50 x 50 lattice 
is shown in Fig. 1. The figure displays our results for ten 
separate recursion lattices, as well as the zero-field result 



-V2- 



for a very large volume lattice, with the eigen- 



mode numbers rescaled by the factor 50^/y, so as to fit in 
the same range on the x-axis as the other ten data sets. It can 
be seen that, at the resolution of this figure, the spectra es- 
sentially all fall on top of one other. The ten separate data 
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Eigenvalue Spectrum, p=18, L=50 
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mode number 



8000 



FIG. 1 : Ten sets of eigenvalue spectra of the operator —D^ — A.o + /n^, 
at (3 = 18, from ten independent 50 x 50 recursion lattices. Also 
plotted, but indistinguishable from the other spectra, is the rescaled 
spectrum of the large- volume zero-field operator — + 

Eigenvalue Spectrum, p=18, L=50 
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FIG. 2: Same as Fig. 1, for the lowest 200 eigenmodes. The closely 
spaced dots are from ten sets of eigenvalue spectra. The "+" symbols 
are taken from the rescaled spectrum of the large- volume zero-field 
operator. 



where we have changed variables, in an axial gauge, from 
gauge field A to field strength B as discussed in the last sec- 
tion. Evaluating in this way the rhs of (38) with f [A,/r[A']] as 
defined by eq. (42), the integration over B is gaussian, and we 
find 



^{R) = -^G{R) 



(52) 



where R = \x — y\ and (no sum over x^y) 



G{R) 
K 



y )yx 



(53) 



Of course, the expectation value of {K )„{K 



-l\ba 



also be evaluated by standard lattice Monte Carlo methods 
based on the D = 3 dimensional Wilson action. A num- 
ber of thermalized lattices are generated by the usual heat 
bath procedure, and K^^ is evaluated on a two-dimensional 
constant-time slice of each three-dimensional lattice. The 
two-dimensional lattices generated in this way will be referred 
to as "MC lattices". They can be thought of as having been 
drawn from a probability weig hting P[U] = "Vl ^p], where 
'I'£,o[f^] is the ground state of the transfer matrix of the D = 3 
dimensional Euclidean lattice gauge theory. 
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sets cannot be distinguished, and the spectrum of — Xq 
looks identical to the (suitably rescaled) spectrum of — at 
large volume. At higher resolution (Fig. 2) some fluctuation 
in the eigenvalue spectrum is observable, and the eigenvalues 
of the lowest-lying modes appear to deviate slightly from the 
zero-field large-volume spectrum. 

Next we turn to the computation of the mass gap. Accord- 
ing to eq. (43), 

(fi> = j DAm2Q[A]P[A] 

« j DAiDA2DA\DA'2 Q[A]P[A,K[A']]P[A'] 

= j DBDA[DA'2 Q[A{B)]P[A{B),K[A']]P[A'] (51) 



FIG. 3: The correlator G{R) computed (i) on two-dimensional lat- 
tice configurations generated from the vacuum wavefunctional by 
the method described in the text; and (ii) on constant-time slices of 
three-dimensional lattice configurations generated by the usual lat- 
tice Monte Carlo method. Lattices generated by the first method are 
denoted "recursion", and by the second as "MC". In each case, the 
lattice extension is 50 sites at (3 = 18. 

Figure 3 shows the data for G{R) at P = 18, averaged from 
a set of ten 50 x 50 recursion lattices, and, for comparison, 
corresponding data averaged from a set of ten 50 x 50 MC 
lattices at p = 18 . Note the very small (~ 0(10^'^)) magni- 
tude of the observable at R — 20, yet even at this magnitude 
there seems to be very little noisiness in the data. Once again, 
this absence of noise is only possible if the variance in the 
K-^K-^ observable is negligible, which supports our original 
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hypothesis. Moreover, the data obtained on recursion and MC 
lattices obviously agree very well with each other 

The mass gap is obtained by fitting the data for G{R) to 
an appropriate functional form, and extracting the exponential 
falloff. Define 



Go{R) = 5"^5^« 
3 



-V2- 



xy 



(54) 



We have seen (Fig. 1) that the spectrum of D~ — Xq is almost 
identical to that of the zero-field Laplacian — V^. With this 
motivation, we introduce the fitting function 



/oW=log 



fl(l 



^MRf 



(55) 



and carry out a two parameter (a and M) best fit of log[G(/?)] 
by fo{R).^ The resulting value for M is an estimate of the 
mass gap. The best fit of the data for G{R) at p = 18 on a 50^ 
lattice by the fitting function exp[/o(/?)] is shown in Fig. 4. 



Correlator fit, |3=18, L=50 
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FIG. 5: Mass gaps extracted from recursion lattices at various lat- 
tice couplings, compared to the 0+ glueball masses in 2-1-1 dimen- 
sions obtained in ref. [16] (denoted "expt") via standard lattice Monte 
Carlo methods. Errorbars are smaller than the symbol sizes. 



The results of extracting M via the best fit of /o to the data, 
for simulations of the vacuum wavefunctional at a variety of 
lattice couplings, are shown in Fig. 5. There we compare our 
values for the mass gap with those reported by Meyer and Te- 
per in ref. [16] (the values for o(p), used in eq. (50), were also 
taken from this reference.) In Table I we list these mass gap 
results, as well as the mass gaps extracted from MC lattices, 
and the "naive" estimate M(0+) — 2m. It can be seen that the 
agreement between the reported values for the mass gap, and 
the masses we have obtained from simulation of our proposed 
wavefunctional (with parameter m fixed to give the observed 
asymptotic string tension), agree within a few (< 6) percent. 
This is a substantial improvement over the "naive" estimate of 
M = 2m, which disagrees with the Monte Carlo results by up 
to 20%. 
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FIG. 4: Best fit (dashed line) of the recursion lattice data for G{R) by 
the analytic form given in eq. (54). 

In an old paper which anticipates the work in this section, 
Samuel [7] argued thatM « 2mo, where mo is the mass param- 
eter in the vacuum state (30) which he proposed. This result is 
obtained if the covariant operator —D^ in (30) is replaced by 
— V^. We believe that a more natural approximation is the re- 
placement of — A-o by — V^, since the lowest eigenvalue in 
the spectrum of each operator begins at zero. Thus the "naive" 
estimate for the mass gap, in our proposal, is M = 2m. 



' The fit.s were carried out by tlie GNUPLOT package, which implements the 
Marquardt-Levenberg fitting algorithm. We have fit the data for log(G(i?)) 
on an L X L lattice in the interval R e [l,L/2]. Errorbars are estimated 
from the variance in mass gaps computed separately, at each p, on ten 
independent lattices. 





mass gap 


(3 




"naive" 


MC 


recursion 


"expt" 






(M = 2m) 


lattices 


lattices 


ref. [16] 


6 


24^ 


1.031 


1.269(5) 


1.174(8) 


1.198(25) 


9 


24^ 


0.627 


0.775(3) 


0.745(5) 


0.765(8) 


12 


32^ 


0.445 


0.562(5) 


0.537(5) 


0.570(11) 


18 


50^ 


0.349 


0.436(3) 


0.402(4) 


0.397(8) 



TABLE I: The mass gaps in D=2-l-l dimensional Yang-Mills theory, 
at a variety of (3 values and lattice sizes L?. Column 3 shows the 
values derived from the estimate M = 2m, and the values extracted 
from G{R) computed on MC lattices are shown in column 4. Column 
5 displays the results extracted from G{R) computed from recursion 
lattices; these are the predictions obtained from numerical simulation 
of the vacuum wavefunctional. All of these values can be compared 
to the mass gaps reported in ref. [16], shown in column 6, which 
were obtained by conventional lattice Monte Carlo methods. 
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VI. VACUUM ENERGY AND CONFINEMENT 

Our proposed vacuum wavefunctional results in a non- 
vanishing asymptotic string tension, via the dimensional re- 
duction argument, for any mass parameter m > 0. In this con- 
text, the question of why pure SU(2) gauge theory confines 
in 2+1 dimensions boils down to why m is non-zero in that 
case, yet m = in the abelian theory. The answer must lie in 
energetics: For some reason the expectation value of (H) is 
lowered, in the non-abelian theory, by having m > 0. 

The calculation of {H) is complicated by functional deriva- 
tives of the kernel K[A]. In this initial study we will simply 
ignore these derivatives, on the grounds that variance of the 
gauge-invariant product K~^K~^ among thermalized config- 
urations has been found, in numerical simulations, to be neg- 
ligible. In fact, this product seems to be remarkably well ap- 
proximated, in any thermalized configuration, by the free-field 
expression Go{R) of eq. (54). This insensitivity to the A-field 
suggests that the variation of K[A] in the neighborhood of ther- 
malized configurations is extremely small, and therefore the 
neglect of functional derivatives of K might be justified. How- 
ever, we are not as yet able to quantify the actual error which 
is made by dropping those derivatives. 

Writing 'Pq = exp(— /?[A]) where 



R 



= -^J d^x<fy B" [xW^B^iy) (56) 



we find 



where 



//^o=(7b-r, + ^| 



(57) 



7b 



5A^ 



T = ^ f^^ dR 

' 2 J 8Al{x) 5Al{x) 



(58) 



Carrying out the indicated functional derivatives of R, but 
dropping terms involving functional derivatives of the kernel 
K leads to 



To 



J d^xd^y8{x-y)i-Dy^K^ 



= iTr 



(59) 



and 



7i = ^1 dhd^yB^x) 
= d^xd'yB^ix) 



-I ab 



Xq — 



B\y) 



ab 



xy 



(60) 




-D- 



^J-lfi-Xo + m^ 

% I d'^xd^yB"(x)( , , 



1 



{Xo — m ) X 

ah 



B\y) 



xy 



Expanding B{x) in eigenstates of —D^ 

B-{x)=Y,bnm 



the second term on the rhs of eq. (61) becomes 

1 

^"V 'X„-Xo + m^ 



2nd term = (^o - m^) -jH^l 



(61) 



(62) 



(63) 



In the previous section, it was found that the eigenvalue spec- 
trum {Xn} is almost unchanged from one equilibrium lattice to 
the next. Then, in the VEV shown in (61), we may replace 

by its VEV with the {Xn} fixed, which is ^g^\/X„ — Xo + m^. 
Then 



2nd term = ^{Xo — rn^)^ 



1 



t \/Xn-Xo + rrP- 




UXo-n?)Tx 



^/-D^-Xo + m^ 



(64) 



iTr- 



Xq - 



^-D^-Xo + nfl ^-D^-Xo + nfl I 



Defining 



■Xfi—XQ 




Xo — 



(65) 



(66) 



(67) 



Suppose that the expectation value of the eigenvalues X„ 
were independent of m^, with zero variance, as in the free 
theory. Setting d{H) /dn? = 0, the minimum vacuum energy 
is obtained trivially, at = (A^). In the abelian free-field 
limit we have A-o = 0, so m = at the minimum and the theory 
is not confining. In the non-abelian theory, in contrast, Xq > 0, 
so = Ao > at the minimum, and confinement is obtained. 
Of course, this simple result neglects both the m^-dependence 
of the eigenvalue spectrum, as well as contributions arising 
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Rayleigh-Ritz Calculation, p=6, L=16 
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FIG. 6: Vacuum energy (H) of eq. (67), per lattice site, computed at 
a variety of mass parameters monal6xl6 lattice at lattice coupling 
p = 6. 



from functional derivatives of the kernel. The situation can be 
improved on somewhat, at least regarding the dependence, 
by a numerical treatment. 

A Monte Carlo evaluation of the energy density (H) /L? as 
a function of m, for P = 6 and L= 16 and (H) as given in eq. 
(67), is shown in Fig. 6. The minimum is away from zero, at 
roughly m = 0.3. This gives a string tension which is a Uttle 
low; the known string tension of the Euclidean theory at P = 
6 would require m = 0.515. This quantitative disagreement 
should not be taken too seriously, because the estimate for 
vacuum energy on which it is based, eq. (67), is of unknown 
accuracy. Once again, in deriving (67), we have neglected 
some terms deriving from functional derivatives of the kernel. 
Even assuming, as we have, that those contributions are quite 
small (and this has not been shown), they could still have a 
large effect on the position of the minimum of a rather flat 
potential. The main point of this section is not to obtain m 
with any degree of accuracy (although that would have been 
desirable), but rather just to see that a non-zero value of m, 
which implies both confinement and a mass gap, is the natural 
outcome of a variational calculation. 



Vn. OTHER PROPOSALS 

There have been other approaches to the Yang-Mills vac- 
uum state in 2+1 dimensions. In particular, the vacuum wave- 
functional proposed by Karabah, Kim, and Nair (KKN) in ref . 
[11] has some strong similarities to ours, and the method we 
have developed for numerical simulation can be applied to the 
KKN vacuum state, as well as to our own proposal. This ap- 
plication is important, because we would like to test the claim 
that a string tension can be derived from the KKN state which 
agrees, to within a few percent, with the continuum limit of 
string tensions extracted from lattice Monte Carlo [17]. 

The KKN approach is formulated in terms of gauge- 
invariant field variables first introduced by Bars [18], and the 



idea is to solve for the ground state of the Hamiltonian, in 
these variables, in powers of the inverse coupling 1/g'^. To 
lowest order, when re-expressed in terms of the usual A-field 
variables, their state has the dimensional reduction form 



l*!,"^ = exp 



where 



2% 



(68) 



(69) 



and Ca is the quadratic Casimir for the SU(N) group in the 
adjoint representation. Because this state has the dimensional 
reduction form, the corresponding string tension is easily de- 
duced. In lattice units, for the SU(2) group, the predicted 
string tension is 



j(0) , 
'KKN 



6 

7tp2 



(70) 



which is in rather close agreement with the lattice Monte 
Carlo results. 

However, the state ^q"' is only the first term in a strong- 
coupling series. As it stands, it implies an infinite glueball 
mass in 2+1 dimensions, and it cannot be even approximately 
correct at short distance scales. The question is whether inclu- 
sion of the higher-order terms in the series, which are neces- 
sary in order to have a non-zero correlation length, will affect 
the long-distance structure, and move the prediction for the 
string tension away from the desired value. KKN resum all 
of the terms in the strong-couphng series which are bilinear 
in their field variables, and when this expression is converted 
back to ordinary A-field variables, their resummed vacuum 
state has the form 



'Po ~ exp 



1 



m 



B\y) 



xy 



(71) 



This state is gauge non-invariant as it stands, and for that rea- 
son must be incomplete. However, KKN argue that the further 
terms in the strong-coupling series, involving higher powers 
of the field variables and their derivatives, supply the extra 
terms required to convert the operator in eq. (71) to a co- 
variant Laplacian. So, according to ref. [11], the vacuum state 
when re-expressed in ordinary variables has the form 



exp 



\/—D'^ + nP-+m 



ab 



B\y) 



xy 



(72) 



In this form, the KKN vacuum state is amenable to the numer- 
ical methods described above. 

At this point we see that there may be trouble ahead for 
the previous string tension prediction. The problem is that the 
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p 




<^KKN 




discrepancy 


9 


24^ 


0.0340(4) 


0.0261(2) 


30% 


12 


32^ 


0.0201(6) 


0.0139(1) 


45% 



TABLE II: A comparison of the string tension Ckkn calculated nu- 
merically from the Karabali-Kim-Nair vacuum wavefunctional (72), 
by methods developed above, with the values of the string tension 
Cmc inD = 3 dimensions, computed by standard lattice Monte Carlo 
methods in ref. [16]. 

coefficient l/(4mg^) in the dimensional reduction form (68) 
is only obtained if the lowest eigenvalue Xq of the covariant 
Laplacian would be zero. We know that this is not the case. 
The effective long-distance wavefunctional, gaussian in the 
fi-field, is obtained as before via a mode cutoff in the B-field, 
and the actual estimate for the KKN string tension, according 
to dimensional reduction, is 

<^KKN = ^ ("^ + Vh)+irfi'^ (73) 

A non-zero Xq will certainly move the predicted string ten- 
sions Okkn away from values given in eq. (70); the question 
is by how much. This can only be determined, at any given p, 
by numerical simulation. 

In Table II we display our results for the string tension 
(^KKN, obtained by evaluating eq. (73) in the vacuum state 
(72) by the methods developed in this paper. It is clear that 
there is a very substantial discrepancy between the predicted 
string tension Okkn and the string tension Gmc, obtained by 
standard Monte Carlo methods in ref. [16]. The disagreement 
becomes disastrous if A,o actually diverges, in physical units, 
in the continuum Umit. In that case the percentage discrep- 
ancy at P — > oo will be infinite. The only way out, that we can 
see, is if eq. (72) is for some reason not the true resummation 
of the KKN strong-coupling expansion. 

An approach which is closely related to that of KKN, rely- 
ing on the same change of field variables, has been followed 
by Leigh, Minic, and Yelnikov (LMY) in ref. [12]. This again 
results in an expression for the vacuum state which is the ex- 
ponential of a biUnear term, with field variables connected by 
a field-dependent kernel. Whereas KKN perform a partial re- 
summation of the strong-coupling series to arrive at their re- 
sult, LMY rely on a conjectured operator identity (eq. (56) of 
ref. [12]) to derive a differential equation for the kernel. The 
hope is that this gives an exact expression for the bilinear term 
in the wavefunction (of course there must be other terms also, 
because the resulting expression for the vacuum is not an ex- 
act eigenstate of the Hamiltonian). Since the derivation reUes 
on a certain conjecture, the justification for the LMY wave- 
functional so far lies in its predictions. 

On the one hand, a glueball mass spectrum resulting from 
the LMY vacuum state has been derived, and this spectrum 
appears to be in very good agreement with existing Monte 
Carlo data. On the other hand, as in the KKN case, the string 
tension (same as (70)) and the spectrum are arrived at by ne- 
glecting the field-dependence of the kernel, which involves 
a holomorphic-covariant Laplacian. We have seen above that 



neglect of the field dependence of the kernel can be dangerous, 
and we think it likely that inclusion of this field-dependence 
will affect the LMY string tension and spectrum significantly. 
It may be possible to use the methods developed here to go be- 
yond the zero-field approximation for the kernel, as we have 
done for the KKN state, to get a better idea of the true predic- 
tions of the LMY state. This is left for future investigation. 



VIII. THE PROBLEM OF N-ALITY 

The Casimir scaling of string tensions is inevitable for the 
lattice Yang-Mills action in two spacetime dimensions, and 
therefore this scaling, out to infinite charged source separa- 
tions, seems to be a consequence of dimensional reduction to 
two dimensions. This feature cannot be true for the asymp- 
totic string tension in 2-h1 and 3h-1 dimensions, except in the 
Nc = °° limit. Asymptotic string tensions in D=2h-1 and 3h-1 
dimensions must depend only on the N-ality of the charged 
source, due to color screening by gluons. The absence of color 
screening in £) = 2 dimensions can be attributed to the fact 
that a gluon has D — 2 physical degrees of freedom in D di- 
mensions. In two dimensions there are no physical degrees of 
freedom corresponding to propagating gluons. If there are no 
gluons there can be no string-breaking via dynamical gluons, 
and hence no N-ality dependence. 

However, the vacuum state of a d + 1 -dimensional gauge 
theory in temporal gauge does not, in general, describe a d- 
dimensional Euclidean Yang-Mills theory, despite the fact that 
each is expressed in terms of a gauge-invariant combination 
of i/-dimensional vector potentials. For example, the vacuum 
state of the 2 -|- 1-dimensionaI abeUan theory, shown in eq. (6), 
describes the ground state of a theory of free, non-interacting 
photon states with a global SU(2) invariance. Our proposed 
vacuum state in eq. (1) interpolates between a theory of non- 
interacting gluons at short distances, and the dimensional re- 
duction form (28) at large scales. If this is the correct vacuum, 
then at intermediate distance scales it describes the ground 
state of strongly interacting gluons with physical degrees of 
freedom; these gluons are free to bind with an external source. 
In that case, the Minkowski-space picture of string-breaking 
via gluon pair production should somehow carry over to N- 
aUty dependence for Wilson loops evaluated in the vacuum 
state at a fixed time.^ 

At present this is only an optimistic speculation, but the fol- 
lowing observation may be relevant: It is possible to compute 
the ground state ^o[f^] in strong-coupling Hamiltonian lat- 
tice gauge theory, and to identify the term in that ground state 
which is responsible for color screening. From the expansion 
of this term in powers of the lattice spacing, we can identify 
the leading correction to dimensional reduction. It turns out 



^ The transition from Casimir scaling to N-ality dependence, due to gluon 

string-breaking effects, is very likely to be associated with a vacuum center 
domain structure, as discussed recently in ref. [19]. Gluon charge screening 
and vacuum center domains are simply two different descriptions, one in 
terms of particles, the other in terms of fields, of the same effect. 
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R[U] = 

Is 





+ Cj 




+ ^2 






+ C3 





+ larger contours 



FIG. 7: The first few terms in the strong-coupling expansion of the 
lattice vacuum state ^o[U], with R[U] = log{^o[U]). 



FIG. 8: How 1x2 rectangles in R[U] screen an adjoint Wilson loop. 
The adjoint Wilson loop (in this case with extension 4x5 lattice 
spacings) is denoted by a heavy solid line. The overlapping 1x2 
rectangles are indicated by (alternately) light solid and light dashed 
lines. The integration over lattice link variables yields a finite result, 
leading to a perimeter-law falloff (eq. (77)) for large adjoint loops. 



that this leading correction has the same form as the leading 
correction to dimensional reduction that is found in the pro- 
posed vacuum state 'Pq [A] ■ 

Denote the lattice vacuum state by M^o[f^] = sxp[R{U)]. 
A strong-coupling technique for calculating R{U) in Hamil- 
tonian lattice gauge theory was developed in ref. [2]. In 
this expansion R{U) is expressed as a sum over spacelike 
Wilson loops and products of loops on the lattice, as indi- 
cated schematically in Fig. 7. The coefficient c, multiply- 
ing a contour constructed from (or filled by) rip plaquettes 
is proportional to (P^)"^. For SU(2) lattice gauge theory in 
D = 2+1 dimensions, the first few coefficients cq,c\,C2,cj, 
of the strong-coupUng series for R[U] were computed in ref. 
[20]. The various terms in can be expanded in a power 
series in the lattice spacing a, and for smoothly varying fields 
it is found that [20] 



Wo[U]=exp 
where 



^ J (fx (aKoB^ - a^K2B(-Z)2)B - 



= 5C0 + 2(ci -t-C2 + C3) 
K2 = tCI 



(74) 



(75) 



and coefficient co is (9(p^), coefficients c\,C2,cj, are O(P^). 
There are several points to note, in connection with eq. (74). 



First, dimensional reduction is associated with the term pro- 
portional to Ko, which receives contributions from all four 
terms shown in Fig. 7, but the leading correction to dimen- 
sional reduction, in the term proportional to K2, comes from 
the 1 X 2 loop in proportional to ci. This is the contour 
which couples B (rather than B^) terms in neighboring plaque- 
ttes. Secondly, it is not hard to see that the 1 x 2 loop in R\U] 
gives rise to color screening. Consider evaluating a spacelike 
Wilson loop in the adjoint representation 

Wadj[C] = joU Tv[Uadj{C)Yi>l[U] (76) 

There is a non-zero contribution to the rhs of eq. (76) which 
comes from lining the perimeter of the adjoint loop with over- 
lapping 1x2 rectangular loops, as shown in Fig. 8, deriving 
from the power series expansion of \U] . For a rectangu- 
lar loop of perimeter P{C) this diagram gives a perimeter-law 
contribution 



ci\m-A 
2) 



(t 



(77) 



to Wadj{C)^ Thus, the same term that gives the leading cor- 
rection to dimensional reduction is also responsible for the 
screening of adjoint loops. Finally, we note that this leading 
correction, proportional to K2, comes in with a negative sign 
relative to the B^ term. 

Now let us consider the leading correction to dimensional 
reduction in the proposed vacuum state ^o[A] of eq. (1). The 
dimensional reduction term was given in eq. (35), and is 
quadratic in B"'"'^. The definition of B*''"" in eq. (34) involves 
a mode cutoff rimax, chosen such that AX = — Xo<^ m^, 
and the first correction to dimensional reduction comes from 
terms in the vacuum wavefunctional of order (A,„ — ?Lo)/m^, 
with n < Umax- These are obtained from the 1/m^ expansion 



1 



1 

m 



2m^ 



(78) 



Taking the second term in the rhs into account, the part of the 
vacuum wavefunctional which is gaussian in 



exp 



m J 



^lowgslow Qslow 



-D^-X 



2rrfl 



gslow _j_ 



(79) 

where the ellipsis indicates higher powers of the covariant 
derivative. We note the similarity of eq. (79) to the strong- 
coupling expression (74). In particular, there is in both cases 
a relative minus sign between the first and second terms. 

The fact that the element responsible for color screening in 
•Polf^] generates, in a lattice spacing expansion, the B{—D^)B 
term coupling B fields in neighbouring plaquettes, is a hint 



^ Generalizing to an SU(N) theory, it is not hard to show (cf. ref. [2]) that 
Co ~ l/s'*^' '^'1 ^ i-/g^N^, and that the perimeter-law contribution shown 
in Fig. 8 is down by an overall factor of relative to the leading area- 
law contribution, as it should be. 
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that it is this term which might be responsible for the color 
screening effect.^ If so, the presence of a very similar correc- 
tion to dimensional reduction, found in 'Pq [A] , would presum- 
ably give rise to the same effect. 

Of course, it is also possible that the vacuum state (1) is 
simply incomplete, and must be supplemented by some addi- 
tional terms which are responsible for color screening. Corn- 
wall [21] has recently conjectured that the dimensional re- 
duction form (28) of the vacuum wavefunctional must be al- 
tered by the addition of a gauge-invariant mass term, imple- 
mented through the introduction of a group- valued auxiliary 
field ^(x), i.e. 



'P[A,<I>] =exp 



- 1 d^'x {ciTr[Fi^]+C2Tr[^-^Dk^]^} 



(80) 

The exponent of this state is stationary around center vortex 
solutions, suggesting a vacuum state dominated at large scales 
by center vortices. This would presumbably solve the N-aUty 
problem. At the moment, however, we lack any direct moti- 
vation from the Schrodinger wavefunctional equation for the 
existence of such a mass term. 

For a discussion of the N-ality problem in the context of the 
KKN approach, see ref. [22]. 

Another type of contribution which is expected to exist in 
the static quark potential is the Liischer — Jt(D — 2)/2AR term. 
We have no insight, at present, as to whether or not this term 
can be generated by the proposed vacuum state of eq. (1). 



IX. CONCLUSIONS 

Our proposal for the ground state of quantized Yang-Mills 
theory, in D = 2 + 1 dimensions, has a number of virtues. 
Apart from agreeing with the ground state of the free theory in 
the appropriate Umit, which is a natural starting point for any 
investigation of this type, we also find agreement in a highly 
non-trivial limit, where the Yang-Mills Schrodinger equation 
is truncated to the zero modes of the gauge field. In addition 
we find, surprisingly, that our vacuum state is amenable to nu- 
merical investigation, despite its very non-local character. 

We believe that this vacuum state may provide some insight 
into the origins of confinement in a non-abelian theory, and 
the precise relationship between the mass gap and the string 
tension. Confinement arises here via dimensional reduction, 
as proposed long ago in ref. [1], and this reduction is obtained 
if the mass parameter m in the vacuum wavefunctional is non- 
zero. We have seen that m ^ is likely to lower the vacuum 
energy, in 2+1 dimensions, and this is related to the fact that 
in a non-abelian gauge theory the lowest eigenvalue Xo of the 
covariant Laplacian is non-zero. The relation between m and 
the asymptotic string tension in 2+1 dimensions is simple, i.e. 



' In fact, apart from an overall sign, the B{—D^)B term looks like the kinetic 
term of a scalar field in the color adjoint representation in two Euclidean 
dimensions. Matter fields of that type can, of course, screen adjoint Wilson 
loops. 



a = 3m/4p, and if the parameter m is chosen to produce the 
string tension known from earlier lattice Monte Carlo studies 
[16], then we find that the mass gap extracted from an ap- 
propriate correlator yields a value within 6% of the mass gap 
obtained by standard lattice Monte Carlo methods. 

The most important unresolved question concerns higher 
representation string tensions. At issue is whether corrections 
to the simple dimensional reduction limit will convert Casimir 
scaling to N-ality dependence, as we have speculated in the 
previous section, or whether some additional terms (such as a 
gauge-invariant mass term [21]) are required. It would also be 
worthwhile to extend our considerations to 3+1 dimensions, 
and to excited-state (glueball and flux-tube) wavefunctionals. 
These possibilities are currently under investigation. 
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APPENDIX A: 3+1 DIMENSIONS 

Although in this article we are mainly interested in the 2+ 1 
dimensional case, it is worth pointing out that the discussion 
in section II can be extended to 3+1 dimensions. Define 



53 = Ai-Ai+A2-A2+A-i-A-i 

C3 = (AiXA2)-(Ai XA2) + (A2XA3)-(A2XA3) 

+(A3xAi)-(A3xAi) 
O3 = [Ai-(A2xA3)f 



The zero-mode Yang-MiUs Hamiltonian is 



H 



1 32 



2 V dAldA^ 



(Al) 



(A2) 



Again we express 'Pq as in eq. (9), and try to solve H^q 
E(^0 to leading order in V. This time, with 



(A3) 



we find 



7b 



y 



-^^+/C3 

dAldAf^ \ 
2x, (^Cl 3D3 



S3 



(A4) 
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In the large volume limit, the ground-state wavefunction 
will only be non-negligible in the region of the "abelian val- 
ley", where the zero-mode components A i,A2, A3 are nearly 
aligned, or anti-aligned, in color space. For definiteness, take 
the large color component (denoted by upper-case A) of the 
color 3-vectors to all lie in the color 3-direction; i.e. 



(A5) 



and lower-case a denotes the small components. With a <C A 
and VRo ~ 0(1) it follows that, in the abelian valley. 





fl} 










Ai = 






aj 








A. 




. ^2 . 




A. 



1 



(A6) 



where A and a denote the magnitudes of the large (color 
3-direction) and transverse field components, respectively. 
Since C3 and D3 are both Oia"^), the non-zero terms contribut- 
ing to To in eq. (14) are at most of order 1 /V^ and can be ne- 
glected. Therefore ^0 — exp[— y/?o], with Rq as given in eq. 
(A3), solves the zero-mode Yang-Mills Schrodinger equation 
to leading order in V, in the abelian valley region away from 
the origin (Aj^ = => ^3 = 0) of field space. 

The generalization of eq. (1) to 3-h1 dimensions is 



^o[A] = exp[-e] 



exp 



ab 



xy 



(A7) 



Again we consider a comer of configuration space in which 
only the non-zero modes make a significant contribution to 
the wavefunctional, and |^Ap Wq. Then 



as before, with 



M''^=535'"'-A^f 



ab 



(A8) 



(A9) 



For a configuration in the abelian valley, with large compo- 
nents in the color 3 direction as shown in eq. (A5), we find 



Q = IgViAiXAjfiM-'l^iAiXAjf 



4 

\gv{MxAjr 



+- 



{AiXAjf 



(AlO) 



Neglecting the overall coupling and volume factors, the rela- 
tive orders of magnitude of each of the three contributions to 
Q are as follows: 



_ {AiXAj)-{AiXAj) 



'Aa"- 



K2 = 
K3 = 



{AiXAj)\AixAjf 



VS3 A 
{AixAj)\AixAj)^ 

m 



m 



Assume that K2,K3 <C Ki. Then we would have 

C3 



(All) 



(A12) 



and I'd [A], evaluated for large zero-mode gauge field config- 
urations, would agree with the ground state solution of the 
zero-mode Yang-Mills Schrodinger equation in D = 3 + 1 di- 
mensions, at least in the neighborhood of the abelian valley. 
But we have already seen that for the solution of the zero- 
mode equation, the magnitude a of the small components is 
related to the magnitude A of the large components according 
to (A6). From this it follows that the assumption K2,3 <C Ki 
in the abelian valley is self-consistent, and justified at large V 
for m^O. 



APPENDIX B: THE SPIRAL GAUGE 

Since ^o[A] in temporal gauge and 2+1 dimen- 
sions is gauge-invariant under two-dimensional gauge- 
transformations, then it is legitimate to carry out a further 
gauge-fixing in the two-dimensional plane when evaluating 
expectation values 

('Polel^o) =IdA Q[A]^l (Bl) 

In particular, with a complete axial gauge fixing, it is possi- 
ble to change variables from A to field-strength B without in- 
troducing any further constraints or field-dependent Jacobian 
factors, i.e. 



DA1DA2 const X DB 



(B2) 



In higher dimensions, as Halpern has shown [3], this change 
of variables would be accompanied by a delta function en- 
forcing the Bianchi constraints, but in two dimensions these 
constraints are absent. 

The simplest approach is to set Ai(ni,n2) = everywhere, 
where (ni,n2) are lattice site coordinates, and invert the dis- 
cretized version of = d\A2 to determine A2 from B. The 
problem with this is that setting A 1 = everywhere on a finite, 
periodic lattice is more than a gauge choice. Gauge transfor- 
mations cannot, in general, set the A field to zero everywhere 
on a closed loop, and lines parallel to the x-axis are closed 
by periodicity. Thus Ai = everywhere is a boundary condi- 
tion, as well as a gauge choice. Although boundary conditions 
should be unimportant at sufficiendy large lattice volumes, we 
would still like to keep such artificial conditions to a mini- 
mum, while retaining the simphcity of inverting B" = d\A2. A 
compromise is what we will call the "spiral gauge", in which 
we set A = (or link variables U = I2) along all links in a 
spiral around the toroidal lattice. An example, on a 10 x 10 
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FIG. 9: The spiral gauge. Link variables on the solid lines are set 
equal to the identity. 

lattice, is shown in Fig. 9. Along the straight sections of the 
spiral, parallel to the x-axis, we have 

A^{ni + l,n2)=B"{ni,n2)+A^{m,n2) (B3) 

For the bent sections, its slightly different. Referring, e.g., to 
the bent section in Fig. 9 starting at ni = 9, = 1. we have 

A1(9,2) = -B''(9,l)-Af(9,l) 
A^(10,2) = B"{9,2)-A1{9,2) (B4) 



Now suppose we start out with setting ^2(1 , 1) =0. Applying 
the above rules all around the spiral we can get all of the non- 
zero A-field variables from the B-field variables, but in order 
to come back to where we started, with A2(l,l)=0, we have 
to require that 

£B>i,n2)=0 (B5) 

ni,«2 

To enforce this condition, we first generate the 5-field without 
constraint, compute the sum 

S" = B^inuni) (B6) 

and then make the readustment 

B"{nun2)^B"{ni,n2)-^ (B7) 

So we have done two things beyond just fixing the gauge. 
First, the A-field has been set to zero on a single closed spi- 
ral around the toroidal lattice. Secondly, by setting in addi- 
tion A2(l, 1) =0, we have imposed a restriction that the B- 
field on the lattice averages to zero in any given configuration. 
These conditions have been imposed for calculational simplic- 
ity; they are not as drastic as setting A 1 = on all links (which 
sets all Polyakov lines in the x-direction equal to unity), and 
ought to be harmless at sufficiently large lattice volumes. 
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